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C \ We investigate the phase structure of an asymmetric fermion superfluid with inter- and intra- 

species pairings. The introduction of the intra-species pairing mechanism in canonical ensemble 
changes significantly the phase diagram and brings in a new state with coexisting inter- and intra- 
species pairings. Different from the case with only inter-species pairing, all the fermion excitations 
are fully gapped in the region with intra-species pairing. 
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Recently, the fermion pairing between different species with mismatched Fermi surfaces prompted great interest 
in both experimental!]], HI, Q and theoretical 0, [Bf, H, 0, @, EJ E3, EL 12, EH studies. In conventional fermion 
superfluid the ground state is well described by the BCS theory, while for asymmetric fermion superfluid the 



phase structure is much more rich and t he p airing mechanism is not yet very clear. Various exotic phases have 
been suggested, such as the Sarma phase fl4| or gapless phase[IH where the superfluid component is breached by 
a normal component in momentum space[ldj|. the Fulde-Ferrel-Larkin-Ovchinnikov(FFLO) phase fl7j] where the 
Cooper pair is momentum dependent, the phase with deformed Fermi surfaces where the two Fermi surfaces cross 
to each other and the total momentum zero Cooper pairs could form at the crossing node[l8|, and the phase 
separation [2fl l2~fl] in coordinate space where the normal and superfluid components are inhomogeneously mixed. 

For many fermion superfluid systems, the fermions from the same species can form Cooper pairs as well. In 
ultracold atom gas like 6 Li and 40 K, there exist pairings between different elements and between different states 
i-^j | of the same element. In color superconducting phase of dense quark matter [22l|, the quarks of different flavors 
can form total spin zero Cooper pairs, and the quarks of the same flavor can be combined into total spin one 
pairs which describe better the cooling rates of quark stars (23|. In neutron starsj^], proton-proton, neutron- 
neutron and neutron-proton pairings are all possible. Recently discovered high temperature superconductivity of 
1271 . l28j can also be well described by an extended two-band BCS theory where the electrons from 



■ the same energy band form Cooper pairs. 

Since the ratio of the intra- to inter-species pairing gap is already 10% for color superconductivity and can even 
reach 50% for nuclear superfluidity and ultracold atom gas, the competition between the two types of pairings 
' should be strong. Considering that the pairing between fermions from the same species happens at the common 
. Fermi surface, the intra-species pairing will become more favored than the inter-species pairing, when the Fermi 
surface mismatch between the two kinds of species is large enough. Therefore, some familiar phenomena in the 
superfluid with only inter-species pairing may be washed out by the introduction of intra-species pairing, and some 
new phases may arise. For instance, the interesting gapless state may disappear, the inhomogeneous FFLO state 
at low temperature 1^, 2^, 33] may be eaten up by the homogeneous intra-species pairing state, and the mixed 



state of BCS superfluid and normal phase [20l |21| may be replaced by a new phase where the inter- and intra- 
i . species pairings coexist. In this paper we propose a general model to investigate the competition between the 
*^ ' inter- and intra-species pairings. 

We consider a fermion fluid containing two species a and b with masses m a and mj and chemical potentials \i a 
and fit,. The system can be described by the Lagrangian density 



where ^(x) are fermion fields with i — a,b and (pseudo-)spin <r =t,4, the coupling constants g,g a and gb 
controlling respectively the inter- and intra-species pairings, are negative to keep the interactions attractive. The 
Lagrangian has the symmetry U a (l) ® Ub(l) with the element U(9 a ,9b) denned as U(9 a ,6b)4>a = e l ° a ip° and 
U(e a ,e b )W =e ie ^. 

We introduce the condensates of a-a, b-b and a-b pairs, <I> Q = —g a (ip^ij}l},§b = — 9b{i ) i' , l ) b) and $ = 
— Server' (V'b V'a)- To incorporate the FFLO state in the study, we assume that the condensates are in 
the form of $ (x) = A Q e 2iq < l ' x , 3>&(x) = A f ,e 2,;qi, ' x and $(x) = Ae i(qa+qi,) ' x with A a , A b and A being independent 
of x. Obviously, the translational symmetry and rotational symmetry in the FFLO state are spontaneously broken. 
Note that, for the sake of simplicity, the FFLO state we considered here is its simplest pattern, namely the single 
plane wave FFLO state or the so-called FF state. 

The key quantity of a thermodynamic system is the partition function Z which can be represented by the path 
integral 

Z = II / ZHW? exp ( J drd^C) . (2) 
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By performing a gauge transformation for the fermion fields, Xa = e _iqa ' x -0^ and Xb = e~ lq6 ' x V>b , the path 
integral over xZ b m the partition function Z of the system at mean field level can be calculated easily, and we 
obtain the mean field thermodynamic potential 
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in the imaginary time formulism of finite temperature field theory where uj n = (2n+ 1)ttT is the fermion frequency 
p is the fermion momentum, and the inverse of the Nambu-Gorkov propagator can be written as 
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with ef = (p ± q l ) 2 / (2mi) - fa. 

The Lagrangian |T|) is non-rcnormalizable, a regularization scheme should be applied. We use the s-wave 
scattering lengthe regularization which relates the bare coupling constants g,g a and gb to the low energy limit of 
the corresponding two-body T-matrices in vacuum by 
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with the s-wave scattering lengthes a a ,ab and a. Such a scheme is reliable in the whole region of interacting 
strength, and therefore it is possible to extend our study to the BCS-BEC crossover, although we focus in this 
paper only on the weak coupling BCS region. For simplicity, we have taken the same particle masses m a = nib = to 
and will take all the condensates as real numbers. 

One should note that if the relative momentum q a — qb is large enough, the uniform superfluid would be unstable 
due to the stratification of the superfluid components characterized by $> a (x) and $b (x) . The case here is analogous 
to the multi-component BEC. To avoid such a dynamic instability 31, 32, 3[|, we choose q a = q& = q 3^] • After 
computing the frequency summation and taking a Bogoliubov-Valatin transformation from particles a and b to 
quasi-particles, the thermodynamic potential can be expressed in terms of the quasi-particles, 
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where are the quasi-particle energies 
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with e±, Se and 6a defined as 



e 2 ± =[(e+e-+A 2 Q + A 2 )±(e+er4 
<ye=(e+-e-)/2 = (e+-e 6 -)/2, 
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A a 7^ and Ab 7^ correspond, respectively, to the spontaneous symmetry breaking patterns U a (l) ® Ub(l) — > 
Ub(l) and J7 (l), and A^O means the breaking pattern U a (l) <8> Ub(l) — » J7 a -b(l) with the element J7(#) defined 
as ?7(0X = e ie C and U{9)ij)% = e^Vb • 

The condensates and the FFLO momentum as functions of temperature and chemical potentials are determined 
by the gap equations, 
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and the ground state of the system is specified by the minimum of the thermodynamic potential, namely by the 
second-order derivatives. 

To see clearly the effect of the mismatch between the two species a and b, we introduce the average chemical 
potential ll = (ii a + /i&)/2 and the chemical potential mismatch <5/i = — /i a )/2 instead of /i a and lib- Without 
loss of generality, we assume lib > We choose ppa,pFa a ,PFab as the free parameters of the model, where 
Pf = \/2mfi is the average Fermi momentum, and set pfcl — —0.58 and pF&a — PFO-b — 0.73 pfcl in the numerical 
calculation. We have checked that in the BCS region of < pf|o|, PiHa a |, PF|a&| < 1 and < \a a \ = \ab\ < \a\ to 
guarantee weak interaction for a-b pairing and more weak coupling for a-a and b-b pairings, there is no qualitative 
change in the obtained phase diagrams. Note that, in the symmetric case with Sit = 0, from a direct integration 
of the gap equation, one obtains the famous result [351] Ao — 8/j,e -2 exp [— 7r/(2pi?|a|)] for the inter-species pairing 
gap Ao at zero temperature. 

We first consider systems of grand canonical ensemble with fixed chemical potentials. The phase diagram in 
T — Sfi plane is shown in FigfT] The left panel is the familiar case without intra-species pairing 36] . Both thermal 
fluctuations and large chemical potential mismatch can break the Cooper pairs. The homogeneous BCS state can 
exist at low temperature and low mismatch, and the inhomogeneous FFLO state survives only in a narrow mismatch 
window. The phase transition from the superfluid to normal phase is of second order, and the transition from the 
homogeneous to inhomogeneous superfluid at zero temperature is of first order and happens at <5/i c = Ao / V% [3|| . 
When the intra-species pairing is included as well, see the right panel of Figfl] the inhomogeneous FFLO state 
of a-b pairing is eaten up by the homogeneous superfluid of a-a and b-b pairings at low temperature, just as we 
expected, and survives only in a small triangle at high temperature. The phase transition from the a-b pairing 
superfluid to the a-a and b-b pairing superfluid is of first order. From the assumption of rib > /i a , the temperature 
to melt the condensate A;, is higher than that to melt A a , which leads to a phase with only b-b pairing. Since 
[Lb increases and /j, a decreases with mismatch Sli, the region with only b-b pairing becomes more and more wide 
when Sfi increases. Note that, for systems with fixed chemical potentials there is no mixed phase of inter- and 
intra-species pairings, and the situation is similar to a three-component fermion system (37^. 
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FIG. 1: The phase diagram of an asymmetric fermion superfluid in grand canonical ensemble. Sfi is the chemical potential 
mismatch between the two species, and the average chemical potential is fixed to be fi — 50Ao with Ao being the corre- 
sponding symmetric gap at S/j, — 0. The left panel is for the familiar case with only inter-species pairing, and in the right 
panel the intra-species pairing is included as well. 



For many physical systems, the fixed quantities are not chemical potentials /i and /ib but particle number 
densities n a = — <90/<9/i a and rib = —dd/dfib, or equivalently speaking the total number density n = rib + n a 
and the number density asymmetry a = 8n/n with Sn = rib — n a . Such systems are normally described by the 
canonical ensemble and the essential quantity is the free energy which is related to the thermodynamic potential 
by a Legendre transformation, J- = Q + fi a n a + Mfc n b = ^ + / in + SliSti. 

It is easy to prove that the gap equations in the canonical ensemble are equivalent to ([9]) in the grand canonical 
ensemble and give the same solutions. While the candidates of the ground state of the system are the same in 
both ensembles, the stability conditions for the two ensembles are very different and lead to much more rich phase 
structure in systems with fixed number densities. Our method to obtain the phase diagram at fixed total number 
density n is as the following. We first calculate the gap equations © and obtain all the possible homogeneous 
phases and inhomogeneous FFLO phase, then compare their free energies to extract the lowest one at fixed T 
and a, and finally investigate the stability of the system against number fluctuations by computing the number 
susceptibility matrix 

dHl Y t R-%l (10) 
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4 



where Yj, are susceptibility vectors with elements (Yi) m = <9 2 f2/ '(dnidx m ) and R is a susceptibility matrix with 
elements R mn = d 2 il/(dx m dx n ) in the order parameter space constructed by x = (A a ,Ab,A,q). The state 
with non positive-definite x (denoted by x < 0) is unstable against number fluctuations and may be a phase 
separation [38|. 

The phase diagram without intra-species pairing can be found in 1^, 2^, 3^] • As a comparison we recalculate it 
and show it in the left panel of FigJ^l By computing the gap equations © and then finding the minimum of the 
free energy J 7 , we find that the superfluid is in homogeneous state at high temperature and FFLO state at low 
temperature. However, the number susceptibility in the region of low temperature and low number asymmetry 
is not positive-definite, the FFLO state in this region is therefore unstable against the number fluctuations, and 
the ground state is probably an inhomogcneous mixture of the BCS superfluid and normal fermion fluid. The 
shadowed region is the gapless superfluid with Sfi > A where the energy gap to excite qusi-particles is zero and 
the system may be sensitive to the thermal and quantum fluctuations. The phase diagram with both inter- and 
intra-species pairings is shown in the right panel of FigO Besides the familiar phase with only inter-species 
pairing (A ^ 0, A a = A& = 0) and the expected phases with only intra-species pairing (A a , Af, ^ 0, A = and 
Ab 7^ 0, A = A a = 0), there appears a new phase where the two kinds of pairings coexist (A,A a ,Ab ^ 0). In 
this new phase the FFLO momentum is zero and the number susceptibility is negative, \ < 0. Therefore, the 
homogeneous superfluid in this region is unstable against number fluctuations, and the ground state is probably a 
inhomogeneous mixture of these three superfluid components. In the familiar phase with only inter-species pairing, 
there remain a stable FFLO region and an unstable FFLO triangle where the number susceptibility is negative 
and the system may be in the state of phase separation. As we expected in the introduction, the gapless state 
appears only in the inter-species pairing superfluid, and in the region with intra-species pairing all the fermions 
are fully gapped. 
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FIG. 2: The phase diagram of a mismatched fermion superfluid in canonical ensemble, a describes the number asymmetry 
between the two species, and the total number n is fixed corresponding to n — 50Ao. The left panel is for the familiar 
case with only inter-species pairing. In the right panel the intra-species pairing is included as well. The shadowed regions 
indicate the gapless superfluid. 



In summary, we have investigated the phase structure of an asymmetric two-species fermion superfluid with 
both inter- and intra-species pairings. Since the attractive interaction for the intra-species pairing is relatively 
weaker, its introduction changes significantly the conventional phase diagram with only inter-species pairing at low 
temperature. For systems with fixed chemical potentials, the inhomogeneous superfluid with inter-species pairing 
at low temperature is replaced by the homogeneous superfluid with intra-species pairing, while for systems with 
fixed species numbers, the two kinds of pairings can coexist at low temperature and low number asymmetry. In any 
region with intra-species pairing, the interesting gapless superfluid is washed out and all fermion excitations are 
fully gapped. To finally determine the exact state of the system in the regions with negative number susceptibility 
(x < 0) in FigJSJ further investigation, especially considering other possible inhomogeneous superfluid phases, is 
needed. 

Since the intra-species pairing is significant at low temperature, our result including both inter- and intra-species 
pairings is expected to influence the characteristics of those low temperature fermion systems. For instance, it will 
change the equation of state of the nuclear superfluid (where neutron-proton and neutron-neutron, proton-proton 
pairings play the roles of inter- and intra-species pairings) or the color superconductor of quark matter (where 
up-down and up-up, down-down quark pairings play the roles of inter- and intra-species pairings in two flavor 
case) in compact stars. On the other hand, systems involving two different alkali atoms can be used to study the 
result we obtained. 
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